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Iru an interconnected power systemj,, if the installed; 
capacity exceeds total demaod plus losses,, more than one 
operating strategy exists to meet the load demand in the 
system and hence-it is importsoat to seek the best one 
that resulits in the minimum} operating cost*.. Economic load 
scheduling of hydro~thermal system is much more complex 
than that of a- purely thermal system,, for two reasons# 

First lyV there is negligible operating cost associated 
with hydro stations*. Secondly, the available water resources 
at hydro stations in a given period are limited# Several 
nonlinear programming techniques have been used in the 
past to solve this problem# But, the basic requirements 
of any algorithm, for implementing optimum operation strategies 
in integrated power systems are speed, efficiency and reliable 
convergence# In view of these factors, successive linear 
programming technique is a powerful candidate # In the 
present work, the application of successive linear progra- 
mming technique is made to solve short range scheduling 
problem of a hydro-thermsl system* Considerable reduction 
in storage and computational efforts, is observed* 
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OHAPOIER- 1 

IirEODDGTION 

With the advent of the high speed digital computers > 
the various problems that arise in the power system > have 
been solved to a very great extent • (Ehe early application 
of the digital computers, however, was restricted^ to the 
problems of analysis such as^ loadfi-flow, short circuit and 
stability studies with greater emphasis on realistic system 
modelling, The economic load scheduling of a large scale 
interconnected power system'^ is also one of the power system 
problems that a computer handles in rigrous manner using 
sophisticated mathematical model for its various components! 
In an interconnected power system with several sources, more 
than one operating strategy is possible to meet the load 
demand in the system and hence it is natural for the operator 

to seek the best operating strategy that results in the 
minimum operating costs, This is achieved by optimising a 
pre-specif ied operating criterion which is the cost of fuel 
^t generating stations, with reference to some control 
variables (plant generations, transformertaps etc.) in the 
system, Q3iis is usually referred to as the economic load 
scheduling, and has been the concern of the power engineers 
for a long time;,- 

A typical power system network consists of a large 
number of thermal and hydro plants connected to the different 
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load centres through lossy transmission networks* Aa impor- 
tant objective in the operation of such a power system is the 
generation and transmission of power to meet the system load; 
demand at each instant of time wilda. the minimum fuel cost# 
fhe study of the optimisation problem is of paramount impor- 
tance in a developing country like India with interconnected 
grid systems having both the hydel as well as the thermal 
plants to meet load demands^ 

Hydro— therm ^ optimal load schedul ing is certainly more 
exciting and complex problem than that of purely theimal system 
optimisation problem, for two reasons ♦= _I'irstly, there is 
negligible operating cost associated with hydro stations and 
hence the hydro generations are not decided by cost considerations 
alone# Secondly, the available water resources at hydro 
stations, in a given period of time, are limited and this 
ppses an extra constraint on hydro generations. In view of 
these factors, the load scheduling of hydro— thermal systems 
bee canes an interval problem, the attempt being to determine 
generations over an interval of time so as to minimise the cost 
of generation in that interval, meeting at the same time, the 
load demand and other constraints • Depending upon the duration 
of the period of optimisation, the problem is referred to aa- 
short range or long range problem# Such a distinction, becomes 
assent ial mainly because of the available information about 
the system* 



The earliest formulation and statement of optimality 
conditions for the economic scheduling of hydro-thermal system 
date back to the pioneering work of Chandler,' Dandeno, Climn 
and Kirchmayer [ 3 ] • This uses Lagrangian multipliers to 
append the equality constraints to the objective function* The 
unconstrained optimisation of the augmented objective function 
results in solving a set of equations together with the equality 
constraints * Dandeno [ 4 ] observed, while working with thesa 
coordination equations on the digital computerj that solutions 
were sometimes obtained that dictate generations beyond the 
plant capa3ity and also negative generation of certain plaatSi 
since the constraints regarding the plant capacities, were not 
incorporated in the initial formulation* 

Wije Per era et* al* [ 5] imposes the upper and lower 
bounds on the operating range of each plant in the form of 
inequality constraints and solved the problem, by using 
Pontryagin’s Maximum Principle followed by Kuhn-dlucker 
conditions. 

Sokkappa [ 27] formulated the hydro-thermal scheduling 
problem as a nonlinear programming problem". Transmission 
losses were considered through the use of loss coefficients* 

Por each subinterval of optimisation, the constraint that is 
most likely to be violated, was picked up and a slack variable 
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was associated with this constraint * The gradient of the 
objective function i*e* cost of the generation,, was evaluated 
and the steepest descent imsthod was used to obtain the solution 
of the problem",' starting from a known initial schedule* Such 
a procedure becomes inefficient when the dimension of the 
problem is large* 

Gopala Rao [11 ] has formulated the short range problem 
as an additively separable nonlinear programming problem and 
used lasdon’s decomposition technique [281, to split the 
problem of higher dimension into sub-problems of smaller 
dimensions to be solved iteratively* The main difficulty in 
this method is the initial choice of the dual variable vector 
which pl^s a dominant role in the computational time required^. 

Dynamic programming technique was used by Singh and 
Aggarwal [26] to solve a short range problem. Haao et.al.[29] have 
used continuous maximum principle of Pontryagin to the economic 
operatiom f or hydro-thermal systems* .Oh[3o] has employed the 
discrete maximum principle for the solution of the same problem* 

In all these methods, the computational requirements for a 
realistic size problem become enormously high. 

OUTLIHE OP THE THESIS 

The central theme of the present work is the application 
of the successive dual linear programming technique to obtain 



5 


solutions to the complex problem of optimal scheduling in 
hydro -thermal power systems including the transmission^ losses • 

A brief outline of the work reported in the different chaptera 
is given below. 

!Ebe second chapter introduces the linear programming 
problem and the simplex method due to George B* Bantzig for 
its solution. The relevance of the dual problem has been 
discussed along with the primal -dual relationship • A' brief 
account of the various nonlinear programming techniques has 
also been given. Nonlinear minimisation via linear approach 
has been discussed along with the algorithm. 

iPhe chapter 3 briefly describes the overview of the 
earlier work done in the field of economic load scheduling 
of the hydro-thermal system. A few of the existing formulations 
have been discussed briefly. 

In chapter 4> the formulation and solution procedure 
for the economic load scheduling problem have been 
presented. The nonlinear equations involved in the formu- 
lation have been suitably linearised around a nominal schedule. 
Successive dual linear programming technique is used to solve 
the resulting linear programming problem. A numerical example 
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is solved to illustrate the solution procedure for a three- 
thermal-one-hydro system* Tranmission losses have been 
considered through the use of loss coefficients* Considerable 

reduction in storage and computation effort is observed*- 

iPjjially chapter 5 concludes with the overall summary 
of the work reported^ in this thesis* A few suggestions or 
possibilities for further work in this firea are also given* 
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OPTBIIZAIIOE TBOHEiaiES 

2 DTTHODUQglOE 

The optimal load scheduling of a hydro-thermal system is 
essentially a mathematical programming problem* It is for 
this reason it becomes imperative to study the optimization 
techniques for solving the mathematical programming problems « 

The mathematical programming techniques referred to 
as optimization techniques are useful in finding the optimal 
of a function known as objective (or goal) function of several 
variables under a prescribed set of constraints ♦ Based on the 
nature of the equations involved, the problem can be classified 
as linear or non?-linear programming problem# A linear progra- 
mming problem (IP or IPP) is one in that the objective function 
and the constraints, both are linear function of the decision 
variables • If any one of the functions among the objective and 
constraint equations happen to be non-linear, the problem is 
referred to as a non-lino ar programming (ELP) problem*. 

2 .2 hlEEAR PRQGRAMMBr& (I.P) 

The linear programming type problem was first recognised 
in the l950s by economists while developing methods for the 
optimal allocation of resources i In 1947, G-eorge B*. 

Dantzig [6] gave the simplex method for the solution of IP*. 
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Since then much progress has been made In the theoretical 
development and also In the practical application of IP. 
Amongst these, the theoretical contribution made by TTubn 

and Tucker had led to the development of the du^ity theory 
in IP [9,10] . 

A standard form of LP enables one to use a standard 
algorithm. All IPs can be cast into the following form. 


Pind the X ?? , 3^), that minimizes* 

The function, P(X) = c. x. (2.1a) 

j=1 

subject to the constraints, 
n 

^i* ... , m. .*(2. 2a) 

j=1 

^ j=1 , ... n ,, (2*3a) 


Ih the matrix form, LP can be, written as. 

Find Z, that minimizes* 

P(Z) = o’^X, (2.11,) 

Subject to 

^ ® .. (2.2b) 

and, where a matrix or vector Inequality Is taken to apply 

to each element of the matrix or vector separately, 

2 > 0 .. (2.3b) 

where A is the nx a matrix. 

*Maxlmlsatlon"of a function is the minimisation of negative 
of the function. 
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In all cases of interest, m < n because of m=n and 
none of the equations are redundant, there is only solution 
to AX=B* If m > n, then there would be m-n redundant equations 
which could be eliminated. With m < n, the system is unde-te . 
mined and, in general possessesan infinite number of solutions 
of which the optimal one is sought* Ihe non— negativity 
constraints X >0 are the only Inequalities present in the 
standard form* 

In order to put an arbitrary problem into the standard 
form, two things are done: 

(a) Inequalities are converted into equalities. If a cons- 
traint appears in the form of a ’’less than" type of inequality 
as 

=k1 ^1 + °k2 ^2 + a-lm ’Si - ‘’k 

it can be converted into the equality form by adding a non- 
negative slack variable follows, 

^k1 “*■ ^‘152 ^2 ^kn ^n+l ” \ 

Similarly, the nonnegative surplus variable is subtracted to 
convert a "greater than typo" of inequality into the equality 
form • 

(b) Variables are transformed, so that, they are all nonnegativ' 
If the variables In the original problem are not required to b 
non -negative, they can be written as the difference between 
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two non-negative variables : 

f »l 

, .. (2.4) 

f It 

where Xj > 0 and x^ > 0 . lEhe new nonnegative variables would 
appear in objective function as well as in constraints* But 
the slack or/and surplus variables would appear only in the 
constraints* 

In connection with IP a nuaber of standard definitions 
are as under* 


(a) feasible solution is one that satisfies Eq* (2*2) and 
inequalities Bq. (2^3) 

(b) A basic feasible solution is a feasible solution with no 

more than m. non-zero x . • In other v7ords, it has at least 

J 

(n-m) X . that are zero * The feasible solution x . asso- 
J 3 

ciatod with the basis vector A^ is known as basic feasible 
solution. 

(c) A nondegencrate basic feasible solution has exactly n 
positive Xj • 

J ■ * 

(d) A optimal basic solution is a basic feasible solution for 
which objective function is optimal* 

The simplex method of Dantzig [6] is a powerful scheme 
to get a basic feasible solution* It makes use of simplex 
algorithm [ 6,7>8 ] Assuming that the problem is feasible 
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and h.as a basic feasible solution readily available, we note 
that if ai^ of cost coefficients corresponding to a non- 
basic variable is negative, the objective function will be 
reduced if we increase the associated variable from zero*^ 
Therefore', to improve the basis, we can bring any such into 
the basis* If more than one is negative, we can choose among 
them, and the criterion most used today is to choose the 't* 
for which 

= min[ c .] .. (2*5) 

j 

This choice selects the variable for which the objective , 
function reduces at the greatest rate* It should be noted 
that some other choice occasionally produce a greater 
reduction in the objective function for this step since, 
because of the constraints, the alternative variable may be 
capable of more change* However, the above criterion is usually 
the best, because it is simplest and also because' experience 
has shown that it is overall efficient* 

If = min[ is non-negative there is no variable 
that can be- increased to reduce the objective function? this 
situation is called the optimality condition* 

We now know which variable to bring into the basis 
(assuming that some < O) and we find, that, we must eliminate 
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the variable for which, 


's 


St 


- min 
^it ^ 


0 


[ b/a3_^] 


• • (2 *6} 


and so, we have the simplex algorithm for proceeding from any 
basic feasible solution to an optimal basic solution if one 
exists. When.^all the in the column corresponding to a 
negative are also negative we can increase x^ indefinitely 

without causing any variable to become negatived This is the 
case when the problem has an unbounded solution. 


The simplex method is a two phase algorithm s phase I 
is to find a basic feasible solution if it is not readily 
available and phase II is to find the optimal solution,- 
out of basic feasible solutions which are finite. The simplex 
method is described in the following steps. 

(i) Phase I of the method ; Arrange the system of Bqs, (2,2), 
so that, all- b^ are p-ositive or zero by changing, 
where necessary, the signs on both sides of any of the 
equationsy 


(ii) Introduce to this system a basic set of artificial 

variables, y.j , y^'i^ , y^^^ where each y^ 0 so 

that, it becomes 

“11 *1 + “12 3^2 + “In Xjj + y, = b, _ ^ 

a2i + 0^22 ^2 ^ ••• + 812x1 ^ ^2 (2 #7) 

I 

j . , . 

V ^ %2 ^ = ’’m 

b >0 
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The objective function of (E ►l^to.aa bci 

written as 

Ci Xi + c„ +..••+ c X + (-f ) =0 * • (2 ♦S) 

I 1 2 2 n n ' v / 

(iii) Define a quantity w as the sum of the artificial variables 
w == y-j + y2 + t y^j^ *• (2.-9) 

«m(l use the simplex algorithm to find x^ 2 0 and 

7±.Z 0 which minimises the Wand satisfy the Bqh»(2*7) 

and (2*8)*. Brom Bqh*. (2 *-7) and (2*9) we get,. 


w = d.| 35j + d 2 X 2 t 

.... t d^x^ + w^ 

i* ( 2 * 10 ) 

where 

.... 

• « 

= - (*1l + ®21 

.... + 

;;.(2.ii) 

for i = l’,2, 

0 H 



and : . . 

“ (b-| + b2 ••• + ** (^t"^2) 

In the algorithm, Eq;* (2 *10), is us-ed together with the 
Eqs* (2*7) and (2*8) for the optimisation*. 

(iv) ’w’ is called the infeasibility form and has the property 
that if , as a result of phase I, minimum of w>0, then no 
feasible solution exists for the original liP problem* 

On the other hand, if minimum of w=0, then the resulting 
solution will act as the starting point for the phase II;* 
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(V) of file method : Apply the simplex algorithm 

to the adjusted system of equations at the end of phase I 
to obtain a solution', if a finite one exists which optimises 
the value of objective function* The flow chart is 
given in the appendix* 

Associated with every hP problem, called the primal, 
there exists another liP problem called its dual* !Che genereil 
primal-^dual relations and the correspondence rules to be 
applied in deriving the general primal "-dual relations are given 
table (2^1) and table (2*2) respectively [8]« 


Table 2 *1 


i^rimal problem 


Ooiresp ending Dual Problem 


n 

Minimise f =3^. c .x 

subject to 
n 


^1 j ^ j ~ ^i* j * • 


jjj - 

Ifeximise v = 3^" y. b. subject to 
i=1 

m ■ • • ■ ■ 

5-^ ^i ^ij'^j * » n|+2,*. m 


1 


IQ 


n 


x.>b., i=m-,+1,m^+2,..m 

where 

where 




x^^O, j=1,2, n^ 
and 

Xj unrestricted in sign, 
j = n-j+l , iij+2, •**, n 


yj_> 0 , i = m^+1, m^ + 2, *. m 
and 

y^ xanrestricted in sign, 

i=:1 , 2 , • * • 131j 
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Table 2 ^ 


Primal quantity 


Corresponding dual quantity 


Objective function ; 

Minimise 
Variable 0 


O^Z 


Objective Punction : 

T 

Maximise Y B 

i"^^ constraint (inequality) 


Variable x^ unrestricted in 

sign 
-th 


.th 


± “ constraint (©quality) 


j constraint^ equality) j^^ variable y^ imrestricted in sign 

constraint, b.(inequa- j"^^ variable y^ > 


^ — 1 variaoj-e y. iO 

j j j 


Coefficient matrix A 



Right hand side vector B 
Cost coefficients C 


T 

coefficient matrix A = 

[ A-j l2» ••• 

right hand side vector C 
cost coefficienIFB B 


Since, in general, an additional constraint requires more 
computational efforts than an additional variable in IP, it is 
evident, that, it is ccMputationally more efficient to solve the 
dual problem, in case the primal problem has more number of 
constraints that its dual. This is one of the advantages of the 
dual problem. 
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2 .5 UQHIiIHEAR PROGRAMMING PROBMI 

A non-linear programming problem may be constrained or 
unconstrained* In most of the cases, the problems are constrained 
one* Some of the most powerful and convenient methods of 
solving constrained optimisation problems involve the transfor- 
mation of the problem into one of unconstrained tjspe* 

(The constraints Involved in NIP may be linear or non- 
linear* For problems having linear constraints, large step 
gradient procedures are available like gradient projection 
method of Rosen [ I 9 ] Por problems with non-linear , constraints, 
a set of necessary conditions to be satisfied by the solution is 
obtained by Kuhn and Pucker [ 20] • Making use of theds condi- 
tions, attempts are made to obtain the solution^ In order to 
point out some of the difficulties in such procedures, a little 
more discussion of these conditions is necessary*- Peschon[ 21] 
has applied a gradient technique to satisfy Kuhn and Pucker 
conditions arising in the reactive power optimisation problem* 

J&. computationally better procedure what is known as the 
sequential unconstrained minimisation techniques (SIMP),, has 
been used by Gopala Rao[11] and Sasson [22] • Briefly their 
approaches are discussed below* 
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The netliod of Fiacco and McCormioli: [24] is based on the 
transformation, 

nx, ^ ^ ^ 

iST j=i 

where f(X);is the objective function, g^(X)<0, i=:1 , •• m are 

Inequality constraints, l.(X)=0, p are equality cons- 

J 

traints and rj^. is the penality parameter ♦ 


A set of decreasing r^^ values are selected and F(X,r^) is mini-, 
mised for each one of them. When r^ is zero (.i>c. sufficiently 
small) minimum of 'F ’ approaches the minimum of ’f ^ * Details can 
be had from Eef • [23, 24] •= 

G-radient methods are widely used for MP problems. ■ 

Many of the commonly used gradient methods start with an initial 
guess, search for a stationary point by making changes in the 
decision variables in proportional to measured partial 'deriva- 
tives to obtain the next guess and so on. The choice of stop 
size is critical and the so-called ’optimum gradient ' [14] is 
the best choice. A few recent techniques of static optimisa- 
tion are also compared In Eef. [16] . Other gradient techniques 
like parallel temgents (PAETAE), conjugate gradient [ 18] . 
Fletcher-Powell technique [1? ], etc. use the gradient and 
some other information to find a new direction.. 
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2*4 LIHEASISAgM APPROACHES TO NOELIHEAR PROBLEMS 

Aru attractive approach to some non-linear programing 
problens is to linearise them in a way that permits the use of 
liP methods. The convergence problem and hi^ solution tine 
make HEP techniques not very useful.. A successive IP formula- 
tion is proposed to correct a potentialEy vulnerable system 

<<w 

so as to bring it to an .optimal normal state... This approach 
has been used by Paranjothi [ 25 ] for optimal power flow problem 
and is also used in the present work for the economic load 
scheduling of hydro-thermal system. Handling equality constra- 

..(2.14) 

..(2J5) 

..( 2 . 16 ) 

let be a solution point which may be feasible or 
infeasible. Expand equations (2 *14)? (2 ..15) and (2^1 6) by 
Taylor series expansion* 


ints in IP is not as difficult as in HEP* 
Consider the original HEP problem? 
Minimise 
f(X), 

Subject to 

S j_ (^ ) 0 > i=1 j . • . m 

and 

i^m = 0 , p 
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fCX) ^ f^X^) + Vf(X^)'^(X-X^) ^PCX) .* (2.17) 
g^(X) ^ g^CX^) + ^gj_(X^)®(X-X^) = G^(X) .. (2.18) 

1^(X) +X71^(X^)®(X-X^) = l.(X) .. (2.19) 

Che corresponding liP problem is as follows j 
Cind X that minimises C(X); 

subject to 

G'^(X^ ^ 0 j i=1 j2j .. in 

and ; 

^j(X) = 0 j j=1 j2, ... p 

Che functions ai*e again linearised about the new point,' 
thus obtained, and the process is continued until the desired 
accuracy in the objective function value is observed. Che 
algorithm can be as under : 

(i) Select a starting point . 

(ii) linearise the functions about this point by Csiylor 
series expansion.- 

(iii) Cormulato corresponding IP problem 

(iv) Solve IP problem so as to obtain a new better point 

(v) Check whether the new point is optimal. If so, stop 
the process. Otherwise, go to step (ii) and repeat 
onwards. 



CHAPTER 5 


BOOHCMIC iOAE SCHEDULIHC OP HYDRO-THEEllAL 

SYSTEMS 

5.1 INTRODHCTIOIT 

The problem of ecLonomic load scheduling of hYdro—thermal. 
systems can be bTOadly classified into l-ong and short range 
problems « The long range problem refers normally to an aaanual 
problem* The inflows^ into the reservoirs show an annual 
csyclicity. Eurthermore'j there may be a seasonal variation in 
power demand on the system. This also shows an annual cyclicity 
The optimisation interval of one year duration is thus a goo® 
choice for long range optimal generation scheduling studies. 

In this case, the solution to the sehed^ing problem consists of 
determination of the sanmounts of water to be drawn from the 
reservoirs fen* hydro-generation in each scheduling subinterval 
and corresponding thermal generations to meet the load demand 
over each subinterval utilising the entire amount of water 
available for power generation during the total interval. 

The short range problem which is the objective of the 
present work,, has an optimisation interval of a day or a week. 
This period is normally subdivided, into hourly or half hourly 
subintervals for scheduling purposes • The so3jition to the 
problem will stipulate the amounts of water to be utilised 
aver each day or week (the scheduling subinterval of the long 
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range problem) • Theni the solution to the short range problem 
consists of an optimal plan for utilisation of this water for 
power generation and the corresponding optimal thermal gene- 
rations are determined considering the load demand on the system 
and the constraints imposed on its operation* 

Several attempts have been made in the past to solve 
the economic load scheduling problem of hydro— thermal systems 
using techniques such as Dynamic programming, Pontryagin's 
TnaT-fTpum principle, KIP techniques etc* In the following section, 
a brief overview of some of them is presented. 

3.2 OVEB.VIBW 01 THE BOOKOMIC IQiD SGHESDIJDIKG PBOBltBM 

Ihe earliest formulation is based up-on the short range 
hydro-thermal economic optimisation problem developed by 
Chandler, etial* [ 3 ] » Ihe main cost of operation of the hydro- 
thermal system is the cost of fuel used in the thermal plants 
as the cost of water in hydro plant is negligible compared to 
the cost of the fuel of the thermal plants. Hence, the objective 
is to determine the generation of individual plants, both 
thermal as well as hydel such that, the cost of total generation 
is the •minitmi'mi a y id at the same timeji total demand and losses 
are continuously met. As it is a short range problem, there 
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will not bo any appreciable chaiagc in the level of water in 
the reservoirs during the optimisation interval Ci»e* effects 
of evaporation and rain fall are ignored) and hence the varia- 
tion in head of the reservoir of the hydel plants will not be 
crucial during the interval* However, specified quantity of 
water must be utilised over the inteival at i"^^ hydel plant*. 
The problem is formulf'.ted as follows : 

Eind the thermal generation S., and the hydel generation 
hj^ which are the function of time, such that the cost function 
as defined by, 

T H 

Cj = I C^(S^(t))i dt .. (3.1) 

is minimum, subject to the equality constraint, 

H M 

y~^s,(t) +y^ h.(t) = + PT-(s,h) (3^) 

i=1 

i^* at each instant of time, the total demand EpCt) and the 
total losses Pj^(S,h} which is the function of plant generations,, 
is equal to the total generation of the plants^, and 
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0} 

I j •' ^5.3) 

in a given interval T, specified amount of water is 
tJa 

utilised at i hydel plant* !Dhe different symbols are as^ 

S.(t) - Power generation of the thermal plant 

J 

hj^(t) = Power generation of the i"^^ hydcl plant 

P = Optimisation interval defining the short 

range problem. 

0.(8. (t))= Fuel cost in Rs/hour of the thermal 
plant which is a quadratic function of 
plant generation S.(t) 

u 

yi(hi(t))= Discharge of water in cubic ft ./hour at the 
i'*'^ hydel plant which is a function of hydel 
ge ne r at ion hj_ ( t ) . 

IT = Humber of thermal plants 

M = Humber of hydel plants 

= load demand at any instant 't • . 

Pj^(S,h) = Transmission losses which are function of 
plant generations* 

The nonlineeuc objective function given by equation (3*1) 
nl.n ng with the equality constraints of equations (3*2) and 
(3^3) are transformed into unconstrained, optimisation with the 
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proper choice of the langrange multipliers X snd The 

following coordination equations are obtained for this optimi- 
sation problem* 

fi_+ X-^=A ..(3 4) 

dS. . 

= X '..(3.5) 

for j = 1 , . , . JT and i = 1, *..M 

The coordination equations (3 *4) and (3 *5)> when solved 

together with equations (3*2) and (3*3), give the solution'p 

for S. 5 (t) and h-(t) over the entire interval* The values Y. 

(known as water value) which effectively convert incremental 

dy- 

water discharge (i.e. as a function of plant generation 

into incremental plant costs, determine the amount of water to 
be utilised at each hydel plant# Therefore, must be so 
chosen that the desired amounts of water are utilised, i*e* 
equation (3*3) is satisfied^ It is observed that is constant 
while X is a function of time . Dandeno [ 4 ] while working with 
these coordination equation on the digital computer observed 
that the direct application of these equations results in solu- 
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tions that sometime s, dictate generations beyond plant capa~ 
cities ard also negative generations because the constraints 
regarding plant capacities are not included in the problem 
formulation. 


(Therefore, the problem is extended by imposing upper 
and lower bounds on the operating range of each plant as 
indicated by the inequalities (3.6) and (5.7) 


c,min 

< 


< 

gmax 

» d=U 

... N 

• • 

and 








, min 
*^1 



< 

j^max 

t i='l > 

... M 

• • 


(3.6) 


(3.7) 


where and h. are the minimum and S. and h^ are the 

maximum limits of operation of the corresponding plant. By 
including the inequalities, Dandeno [ 4] doubted about the 
constancy of The preliminary step in the extension is 

the application of Pontrvagin’s maximum principle by Wije 
Perera [5] . By applying this maximum principle followed by 
Kuhn—Tucker conditions, the problem' reduces to maximising the 
Hamiltonian H at each instant of time i.e. 




Maximise H = 0^(S^(t)) - ^ T yj_(hj^('fc)) •• (3*8) 

i=1 
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The maximisation of H at each instant of time is an auxiliary 
problem of static type • The value of should bo so chosen 
that the equation (3*3) is satisfied Y^ actually corresponds 
to the price of water* An increase in Y^ would result in 

•f'Vi 

the lesser water usage at the i hydel plant and vice— versa*. 
Maximisation of is equivalent to the minimisation of the 
function of equation (3*8) as shown below : 


Yi •* (3*9) 


The minimisation of the function as defined by 
equation (3*9) corresponds to the minimisation of the total 
operating cost of the thermal as well as hydel plants » The 
first term in equation (3*9) represents operating cost of the 
thermal plant and the second term, that of the hydel plant 
which is almost negligible in practice , 

Gopala Eao [11] formulated the economic load scheduling 
problem of hydro— thermal system ae sen additively separable 
nonlinear programming problem and used Lasdon’s decomposition 
technique [ 2^ for its sciinticai*.' The optimisation interval is 
subdivided into ’K’ equal intervals such that the load' demand 


Minimise = Maximise (-(H)) 


IT 


M 


3=1 1=1 
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at all the stations remains constant*, lo account for losses, 
A.G. power flov? equations are used* Bus voltage nagnitudes 
and angles constitute the problem variables » Ma-sri mnin and rndni- 
mum limits on active and reactive power generations and maxi- 
mum limits on bus voltage magnitudes are imposed*. 

Iiet theTe be *M‘ hydro and thermal stations and 
be the. total buses*. Then,; in any k interval-, the injection 
equations can be written as 

m 

qC ^3 ^13 ^ •• 

j=1 

m 

Jj_( E^) = ®i ® j ^ij ®ij ~ ^ •• C5-11) 

A— i* i.i — » V 

where E^ ard ^ ^ are bus voltage magniti;de and angle in 
interval 


^ij is the ij^^ element of Y bus, 

*feii 

injections into the system at i bus*. 


and and are 


Then power generation equations will be 


C 6^,B^) = I. ( 6^,1^) + pj. 

( 6^,E^) = J^( 6^,b^) + 

fori = 1, .•.,(M+ir). 
and k = 1, ••», K 


.. ( 3 . 12 ) 
.. ( 3 . 13 ) 
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where and are active and reactive power generations 
"tjii ic <i*vk 

at i bus, and and are active and reactive loads at 

ifcii 

the same bus in the k interval 


Similarly, the load demand eq.uations will be, 

B^) + = 0 .. (3 JH) 

E*") + = 0 .. (3.15) 


for i = M+h+l , m, 
and k = 1 , K 

fhe objective function, i.e.. the cost of generation is 
given as , 


K M+U 

f( E^) = y~* 0 (E^( 6^, e'^)) .. (3,16) 

/ \ Z- A -L -i- 

k=1 i=M+1 


Now, the optimisation problem can be mathematically 
stated as, 

Minimise f ( 

Subject to (i) load demand eq.uations, 

I^( E^) + pJi = a ff (3*17.) 

J^( 6^f E^) + = 0 f (3*18) 

for i = M+N+1 , • • •, m 
and k = 1, •••K 
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(ii) engineering limits on the system 


pinin 

S. P^( 6 

vl 

m 

pmax 

i 

:: (3 '.'19) 


<. 6 

S^] 1 


. (3'^0) 


<- 

j,max 

. i 


.. (3.21) 


for 

i = 1, 

. * • M+N 



and 

k = 1, 

... K 



and (iii) coupling equations at hydro plants 
K 

y " antl’iC 6^1 E^)) = \ *.(3,22) 

k=1 

for i = 1 , , » • M 

where is the volume of water available for power generation 

fhe objective function can be rewritten as, 

K 

f(6^, =\ f'( 6^, E^) .. (3.23) 

I ^ 

k=1 


where 


M+U 


= y 


C.(P?( E*")) 


i=M+l 


.. (3,24) 
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Associating dual variables u>j , •«. with coupling 
equations, the lagrangian function can be written as,' 

K M K 

k=1 i=1 k=1 


• • ( 5 ^ 5 ) 


Now the dual problem will be, 
Maximise h(u) over all u 


where h(u) = mini (6^, E^, u) .. (3^6> 

for k = 1, K 

Ihe Set (known as constraining set) to which the k"^^ sub~ 
problem variables ^ and are to he constrained, is defined 
by the equations ( 3 * 17 ) to (3 •Si), 

Using decomposition technique, the lagrangian function 
can be minimised by solving the following ’K’ subproblems 

for k = 1 , , ♦ • K 

M 

min 

i=1 
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subject to (i) the load demand equations of (3fl7) 

( 3 - 18 ) 

and (ii) the engineering limits of equation ( 3 * 19 ) 
to (3 •21) ,: on the system* 

fo sum:: up, an additively separable nonlinear programming 
problem is split into subprobloms and solved iteratiTely- The 
choice of dual variables is crucial and plays a dominant role 
in the computation efforts required*. 

Singh and Aggaiviral [26] considered a short range problem 
and attempted to solve it through the Dynamic programming 
technique* It assumes, that the thermal plants are run as 
base load plant and therefore, depending upon the demand, 
all the thermal plants will be put into operation one by one 
and afterwards when the demand exceeds, hydel plants that are run 
as peak load plants to take core flf fluctuating loads will be 
put into operation* It does not include transmission losses 
and they are computed for each interval after the schedule has 
been obtained for a demand corresponding to that interval.^ 

This is a serious drawback indeed* furthermore, the Dynamic 
programming becomes Inefficient when the dimension of the 
problem is large ^ 



CHAPTER 4 


EGOMIG LOAD SGHEDUIiIHG OE HTDRO-THEEMAIi SYSTEMS UGLUDIHG 

TEANSMISSIOl LOSSES 

4.1 IHTROLGCTIOH 

The crux of the load scheduling problem is to minimise 
the total generation cost of the power system while satisfying 
all the constraints imposed on the system . Several nonlinear 
programming techniques have been proposed in the literature 
to solve such problems [ 3,26^27,29^30] . While they differ 

in the type of constraints that they can handle and techniques 
used, the general approach is to solve the problem so as to got 
new set of plant generations that reduces the cost of genera- 
tion and repeat the procedure until no further reduction in 
the objective function i*e* the cost of generation is possible# 

In the present chapter, the formulation and algorithm 
to solve the short range load scheduling problem of hydro- 
thermal systems are presented* In a combined hydro— thermal 
system., there is no (almost negligible) cost of generation 
associated with hydro stations* Ho7/cver, there is always a 
constraint on total energy to be generated at hydro stations* 

The total energy available at each hydro station over any given 
period of time is determined by the amount of water available for 
generation* This is why the formulation of hydro— thermal 
problem can not be an instantaneoirs problem as one finds the 
case with purely thermal systems because it is assumed that 
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the fuel supply to thermal systems is continusously available*. 

The short; range problem attempts to allocate the resources 
obtained from the solution of long range problem over a day 
so as to minimise the cost of generation during that day*- 
Iiong— range problem is not a part of the work reported in this 
thesis* Hence, in the numerical example, certain amount of 
water is assumed for utilisation over a day* It is expected 
that one realises more savings if the transmission lino losses 
exe analysed accurately*. Transmission losses have been taken 
into account by loss coefficients developed by Kirchmayer [l ,2 ] 
The uncertainties in load demand during a day are usually negli- 
gible • Hence, it is reasonable to treat load dennnds to be 
deterministic* However, the consideration of head variationi 
may not be crucial for a short range problem* With this 
assimiption, by discratizing the time into reasonable hourly 
or half hourly intervals,, the short-range can be formulated and 
solved by the successive linear programming technique after 
suitably linearising the functions involved about the nominal 
solution pointy 

4 *2 gQEMHEiATIOH 01 THE PROBitBM 

As mentioned earlier, the interval of optimisation 
(usually a day) is discretized into ’K* equal intervals, such 
that in each interval the load at all stations remain constant*. 
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It is ftirther assumed that the solution to the long range 
problem is available, which specifies the hydro energy to 
be generated at each hydro station. The system that we are going 
to study is predominantly a thermal system with one or two 
hydel plants. IChis is because of our basic assumption that the 
thermal plants are operated as base load plants and the hydel 
plants are operated as peak load plants . However, if the 
number of hydel plants are more than one, they can always 
be combined into a single equivalent pliant on the assumption 
that each of the hydro*^stations ^e located at about the same 
distance from the load centred This is because the incremental 
transmission loss is about the same . Thus'^' after obtaining 
the generation schedule of the equivalent hydel plant, we can 
assign this generation to tho various hydel plants in proportion 
to their ratings. 

I»et there be 'H' thermal plants and hydel plant s> 

let be the generation of the thermal plant and hji bo 

the generation of the i^^ hydel plant. Kie objective function 

in the present work very closely follows the formulation due to 

Singh and iggarwal [ 26] . let the cost of generation in any 

interval for a power (thermal) S^ generated at j thermal 

plant be G.(S.)*. If the intervals are not of equal duration, 

J. iJ 

k 

of courSG, tho function 0^ will bo doscribod by different 

tj 

coefficients in each interval, otherwise they do not depend om 
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the interval* Then the overall cost of generation is given by 
K H M 

C^^CSyh, V ) ■^1^1 (iif) .*• (4 -1) 

k=1 i=1 i=1 

In the Bq* (4*1)j Cu is a function of thermal generation 
S^, hydel generation h^ and a fictitious cost of water T^. 

While optimising the Bq* (4*1) under the set of constraints, 

Yj_ is made pre-assigned parameter i*e* some appropriate value 
is assigned to Tj • The significance of Y . is to associate a 
fictitious cost with the hydel generation, otherwise,, optimisation 
pxocess would try to assign all the lo^s to the hydel plants^ 
only * Hence, Yj_ controls the generation of the hydel plants » 
which in turn controls the discharge of the water. The dimen- 
sion of YjL is that of water price i*e* R/cubic ft* 

yfChf) is the discharge of water at i"^^ hydel plant in ci*ic ft 

•bh 

hour for a hydel generation of hj|^ at the i hydel plant* Sic 
higher the value of Y^'^' the lesser the. discharge of water* 

However,, if we do not put any constraint in the optimi- 
sation of the objective function as defined by the equation 
(4*1),. wo would get the trivial solution i*e * zero cost with 
no generation* Therefore, the following constraints must be 

imposed* 
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N 


M 


a=i 


* 2 


1=1 


>^1 

r**^ 3« 




• • (4 *2 ) 


for k=1 , . . . K 


The eqixation (4*2) states that for any interval of time, 
the total generation duo to hydro as well as thermal generation 
must be met by the load demand plus losses in that interval.^ 

The losses are a function of plant generations* They can be 
suitably represented by the loss ooefficient matris. known as ’B' 
matrix [1,2]* 


IT 




M M 






j=1 i=1 


+ 

^ 1 


_ ^ ^ 

/ ^i hj 


j=1 i=1 


^ 3 


M M 



3=1 i=1 





. . (4 *3) 


where b^p^ , b t^ and b^ are the elements of the 
13 ij ij 

submatrices Bjtfpj and B^ such that the matrix or the loss 
coefficient matrix is as under, 
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B = 


B 




B 


!EH 


B 


T 

m 


B. 


HH 


.. (4^) 


Care should be taken that the solution obtained should 
not dictate generation of any plant beyond the plant capacity 
or negative generation’* This is why the following inequality 
constraints miist also be incorporated 


^ . TniTi _k oiaax 
S. < S. < S*7 

3 ” 0 3, 


.. C4*5) 


, , min ^ , k ^ .max 

and h^ ^ hj^^ ^ ^i 


• • (4.6) 


for i=1 f • • • M 
j =1 , • . • H 
ajcid k=1 f • • • K 


However specified amount of water le discharged at 

the i^^ hydel plant over a day. This forms another constraint 

’fctl 

on the hydel generation h^ at the i hydel plant •• 

K 

k=1 



38 


In the following section, it will be shown that the 
problem as formulated here, can be solved by using successive 
linear programming technique which has the advantages of simp- 
licity in formulation, fast in execution and at the same time » 
ftssuraaoe of a solution If one exists. 


4.3 IIEBAR PROGRAMMING PROSPEM POBMUIiAIIQU 

The minimisation of the system generation cost given by 
the equation (4*1) is considered as an objective function* 

The cost characteristics of the different theimal plants and 
discharge characteristics of the hydel plants as a function of 
the pliant generations are known* These are non-linear functions 
and usually represented by a quadratic as shown below, 

Oj(S^) = a. + bj •• (4.8) 

yi(h^) = Pi + •• (4.9) 


where G.(S't) 
j J 


= aost function of j^^ thermal plant in Rs//3R 

'tih. 

= thermal generation of j thermal plant in W 

•f'Vi 

= discharge function of i^ hydel plant in C-H.ft./EIR 
= hydel generation of i^^ hydel plant in Mff 




Pi>. (^i Slid r^ = constants. 


and 
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!Big constants in the equations (4*8), and (4 *9) afire 
assumed to be known. The objective function of the equation 
(4,1) can be rewritten as, 


K N 

k=1 j=1 



M 



Vi.(Pj. + li >"1 + 


.. (4 *. 10 ) 


Our approach would bo here to minimise the cost of 

generation in each interval of time over the total optimisation 

period. While considering the cost associated: with a particular 

k • k 

interval ‘k’, we may drop the subscript 'k' from 3^ and h^ and 
write our objective function for that particular Interval as 
follows 

U M 


f ^ = 






i=1 


i'rom the equation (4«H)> it is evident that 0^ in 
equation (4.l0) is nothing but the sum of ’f^‘ for all K intervals 
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over the total optimisation period*. We logically develop a 
feel that the minimum of is the sum of minimum of ’f^^ for 
all K intervals over the total optimisation period* 

Linearising the scalar cost function of equation (4 ^ ) 

by Taylor series expansion around the initial operating state 
(Sji h^)j, the incremental cost function is obtained as, 

N M 

A8j +\ AH^ .. (4.12) 

4=1 


i. # 



jL 

i=:1 


dfAPi 


r, (4*13) 


where d^ is the incremental cost coefficient and 
is the incremental change in the generation of the plants* 
and A are given as under 

d. = b.+2c. •• ( 4 * 14 ) 

J J J J 


for 3=1 , ... N 

%+l = Yl (% + 2 =^1 1^?) . 

for i=l , . • • M 


.. (4.15) 
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AP^ •• C4..'16) 

for j=1, •'.. ir 

Ap^^^ = Ah^ :: (4;i7) 

for i=:1 , •*, M 

Ihe constraints of the equations (4*5)» (4*6) and (4*2) 
OEsn he rewritten in the following form, after the dropping the 
subscript ’k* 


g. = s.-st 

®3. 1 i 


,min > 


0 




for i=1 • ST 


, , min - « 

S2]J+i = ‘‘i - ‘'l ^ ° 


^ , max^ V, s n 

S2M+i ~ \ ~ ^i - ° 


for i=1 , M 


... (4.^1 8a) 
.. (4J8b) 

C4**4^a) 
;l (4.l9h) 




M 


^2F+2M+1 


h,. - 


3 ^z- - 

d=1 i=1 



P^-Pj.CS^Vlii)=0 (4.20) 


’jm-AM’i 


lift. 
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linGarising the equations (4»18), (4 *19) and (4*20) 
around the ciu?rent operating state by Taylor series expansion, ^ 
we get the new linearised constraints as follows* 


gi = AS^> 3^ - s° 


(4 .21 a) 


For i=1 , 11 


%+i 


^2N+i 


~ 


-AS. > S° - 
1 “* 1 1 


for i=1 , • F 

-hj 

for i=1 , ••• M 
-Ah^ > h? - 

for i=l , M 


.. (4.21b) 


.. (4.22a) 


.• (4.22b) 


F F 


M 


^2F+22fi41 “ 




i=i 

M F 


m=1 


i=l n=1 


M 

-.V: 

"nl 1 

Ill=l 

F M 




P^.P^(S- h^)-^s- -^h! 

. 1=1 1:^1 


(4 '.25) 
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Equation (4.13) togethor with equations (4.21),; (4.22) 
and (4.23) forms a liP problem in incremental variables A S. 
andAlij_* “lain requirement of the IiP technique is that 

the problem variables are restricted to be non~negative » But 
in the formulation proposed here, the incremental variables are 
unrestricted in sign. To overcome this difficulty Dantzig 
presents two methods* In one method,; the unrestricted variables 
are expressed as the difference of two non— negative variables 

This doubles the number of the variables in the problem and 
loads to increased storage requirement to solve IP problem. 

Another method is to shift the unrestricted variables, if the 
minimum or maximum limit of the variables are known, to ensure 
nonHaegativity.. This reduces the number of constraints because 
the constraints on minimum or maximum limits for the problem 
variables are not required. The second method has been 
preferred in this case because, (i) the limits on the problem 
variables are known (ii) storage requirement is reduced and 
(iii) number of constraints gets reduced, Therefore, we define 
a? non-negative vector 2 of dimension (H+M) with elements 
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as follows, 


.,=As. -dsf> 


for j=1, ... N 




for 1=1 , . M 


(4 •24 a) 


.. (4?24b) 


where 


t) «) 

ard Ah;^=h“^ 


• • (4.25 a) 

• • (4.25b) 


With the variables defined by the equation (4.24)j the 
equations (4^1 a) and (4.22a) are redundant » Substitution of 
equation (4.25) in the equations (4.2lb),, (4.22b) and (4^3) yields'jij 


-z. > 


• • (4 .26a) 


for j=1, .. • N 


^ , mm .max 

- ‘'i - **1 


.. (4.26b) 


for i=1 , • • . M 
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N 


M 

'I o 


M ir 




i-1 11=^ 


■ja 


m. 


1^1 


1-2^s« bjH^ 

i=l n=1 



M 


F K? 


-2 > \ ^HH 

. MX. 


Pj,-pj,(s°,h°)-y(i=4\ 


m=1 


L — i 
3:=1 n=1 


JJI 


M N 


M 




(1-2 \ S° b— -2\ ft® 


'n “TH, 


-2) iC bj2j )hl 

•ni / ^ 


,inin 
. « 

1 


in=l 


i=1 n=1 


ia=l 


.. (4.27) 


Substitution of equation (4.25) in the expression for 


incremental objective function given by equation (4.13) yields 

j. 


Af^ = <f d^ 


ginin 




.. (4.28) 


In the equation (4.28), the second term is constant and 
does not contribute to the optimisation process. Iho objective 
function reduces to 


A f ^ Z 


. . (4.29) 
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!Dhen liP probl©ia can be stated as • Determine 2- tiiat 
minimises the objective function of equation (2,29) subject 
to the constraints of equations (4 •26) and (4 •27). 

In the present formulation,, the number of constraints 
exceeds the number of variables . Hence dual IiP problem is 
obtained of this primal problem, to attempt the solution 
because less effort is needed to solve the problem when number 
of constraints are less than the number of variables. Tables 
(2.1) and (2‘.2) in the section (2 .2) of Chapter 2 can be used 
to construct the dual problem 

4.4 SOLUTION PROCEDURE 

We shall now describe the solution procedure using dual 
IiP technique for the economic Ictad schedul i ng problem * 

Step 1 * Assume for all i*s. Its choice is crucial. 
Usually WQ start with low value such as 0^1 • 

Step 2 : Assume an initial feasible value of S® and h^. That 
is, select the initial and h^ such that it satis- 
fies the condition (4^5) and (4*6) 

Step 3 * Set k=1 



47 


Step 4 : Set Pp = Pp 

Step 5 i Obtain B and C vectors for primal problem of equation 

(4 *29) by using the R»H*S* of equation (4*26) and (4»27) 
and equations (4^14) and (4 *15) respectively^ Formu- 
late the dual IiP problem using relations and rules 
given in the Tables (2,.1) and (2*2) of section (2^) 
of chapter 2 * 


Stop 6 
Step 7 


Step 8 
Stop 9 


t Calculate the system cost of generation with S° and h? 

0 i 

: Solve the dual IP problem to get corrections i.e.AS. 

3 

and A h^* 


s I^dato the and h^ as = S? •♦■AS^ and =h?+Ah^ 


Calculate the system cost of generation with S^®^ and 

^new 
^i • 


Step 10» If the difference in the system cost of generation is 

less than the specified tolerance, store S^^ and h^?®^ 
k k j 1 

as and h^ and go to step 11 • Otherwise set 

and h^ = h^®^ and go to step 5 and repeat onwards* 

Step 11: If k is less than K, set k=k+1 S®=S^ and and 

go to step 4 and repeat onwards* Otherwise, go to step 
12 * 
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Step 12 s Check for the discharge of water by using equation 
(4 ,7)., If the discharge is within tolerable band, 

Stop the procedure.* Otherwise adjust and go to 
step 2 and. repeat onwards V 

!Ehe efficiency of the proposed method greatly depends 
upon the choice of the water values i^* As already stated, 

this method is very much powerful for a system which is pre- 
dominantly a thermal one with one or two hydro-stations* 

However, they can be combined into one equivalent source for the 
scheduling purpose, if they are large in number'*- 

4.5 HIMERICAIi MAMELE 

The formulation and the solution procedure discussed in 
the preceding sections is applied to a short range lofd schedu- 
ling problem of_24 hours duration as the total optimisation 
interval.. Ihc total optimisation interval has been subdivided 
into equal 24 intervals of one hour each. (The demand at the 
generating station is assumed to bo constant during each hourly 
interval and at the end of each interval, the demand increases 
or decreases in jumps (or remains constant^ • The system consists 
of four generating plants in all, out of them^ three are thermal 
and one, hydel plant*. !Mio cost discharge characteristics are 
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giron in Table 4 ►I • 


Table 4 .1 


Oost/discharge oh-aracteristics 


lower bound 
in W 


Upper bound 
in MW 


Cost characteristic 


Quantity of water 
to be utilised 
over 24 hour in 
cubic ft • 


50 

60 

50 


200 

170 

215 


(a) OKiermal Plant 
0 ^= 100 + 0 . 13 ^ + 0 . 01 3 ^ 

02 = 1 20+0 .1 32+0 .023^ Rs /hr 
0^=1 50+0 .23^+0 .01 s| 


15 


65 


(b) Hydol Plant 
=1 40+20hj_+0 .06h^ 
Cubic ft ./hr 


18,000 


The loss coefficients are taken to be as under. 


0.0005 

0.00005 

0.0002 

0.00003 

0.00005 

0.00004 

0 .0001 8 

-0 .0001 1 

0.0002 

0 #0001 8 

0 .0005 

-0.00012 

0.00003 

-0 .0001 1 

-0.00012 

0.00023 
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The above problem was solved on.DEG-10 Digital Computer 
at IIT Kanpur» using the formulation and the solution procedure 
discussed in the preceding sections for the load demands given, 
in the Table (4 *2) for each intervals The results of generation- 
schedule of all the thermal and hydro plants were obtained as 
shown in Table (4*2)« Transmission losses were taken into 
account by the loss coefficients* The problem converged for a' 
water value of 0*22. 

4-6 DISCUSSION 

The general formulation proposed in this chapter has 
shown that the economic load scheduling of Hydro— thermal system 
can be done by using successive dual DP approach* The 
equality constraints which arc difficult to handle in NDP 
problems, can be easily handled by this approach*. 

Since the minimum system generation cost for a parti- 
cular sub interval is obtained practically in 1 to 2 iterations, 
this approach is really well suited for large size problems*. 
Number of constraints are suffici^tly less than the actual 
constraints one has to tackle with NDP techniques because 
no constraints are required for minimum (or maximum) limits 
for this type of formulation* This results in saving of 
storage, and computational efforts* For a very bad starting 
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Table 4 .2 


Resul-ts of Generation Schedule 


Inter- 

val(k) 

Demand 

cpj) 

±Il(Kw) 

Optimal 

cost 

(f^) in 
Rs* 

Generation Schedule 

Thennal Plant . 

k ^k k 

Si S 2 S3 

in (Mw) 

Hydro- 

plant 

^1 

Discharge 
water 
(in cubic 
ft) 

Trans-r- 

mission 

losses 

1 

175 

591 49 

5.5 .2 

60 4 

50 4 

154 

453 .'5 

5-2 

2 

190 

612.82 

61 .2 

60 .0 

60.5 

15 .0 

453 .5 

6 *7 

3 

220 

665 44 

73 .6 

60.0 

81 .8 

15.0 

453 4 

104 

4 

280 

810.00 

1 074 

60 4 

117.5 

1 5 4 

453 4 

2p.7 

3 

320 

932*19 

1 31 .8 

79 4 

120*3 

15*0 

453 4 

26 4 

6 

360 

1079 45 

157*1 

81. 8 

142.7 

154 

453 4 

36*7 

7 

390 

1203*80. 

159.3 

96.4 

1.62 4 

15 .0 

453 4 

45.1 

8 

41 0 

1292 .86 

171 .0 

107*3 

162 4 

16.0 

474.5 

46*7 

9 

440 

1421 .63 

189 '.2 

109 4 

164.5 

28*9 

767 4 

51 .6 

10 

475 

1603 .60 

189.2 

134*7 

176 4 

31 *3 

824*5 

57*1 

11 

525 

1 842 *16 

198*9 

1 34 .7 

2024 

55 4 

1436.2 

66 4 

12 

550 

1929.23 

2004 I 49 .O 

204.5 

65 .0 

1693.5 

684 

13 

565 

2080*14 

2004 157.8 

215.0 

65 4 

1693 4 

72.8 

14 

540 

1914*92 

2004 1394 

203.0 

65 4 

16934 

67.2 

15 

500 

1 703 .63 

189*3 

134.7 

180*5 

52 .8 

1363*1 

57*3 

16 

450 

1471 .61 

189*2 

134*7 

•176*9 

51 .3 

824 .6 

57*1 

17 

425 

1347 .01 

171 4 

107*3 

164 .3 

284 

7674 

46-.7 

18 

400 

124844 

1714 

97*5 

162.4 

15.0 

453.5 

46 .0 

19 

375 

1138*74 

157*1 

96 4 

144*7 

154 

453 4 

38.2 

20 

340 

1003 46 

134.6 

79 4 

142*7 

15 4 

453.5 

31 *7 

21 

300 

870 43 

1104 

79 4 

1174 

15.0 

453*5 

22*1 

22 

250 

730 *1 7 

1 04 4 

60 4 

85.9 

154 

453.5 

14.9 

23 

200 

628*61 

71 *1 

60.0 

61 •7_ 

15.0 

453*5 

7.8 

24 

180 

598';3l 

60*1 

604 

504 

154 

4534 

5.6 - 


Total Syst^ cost of generation = Rs»28770»36 
Total discharge of water = 17830V7 Cubic ft* 
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with as zero-, the computGr time required to obtain the 
solution to the numerical problem in section (4»5) is 30 seconds 
OiiI>BG~l0* But for a judicious choice of Yj_ as 0«2, it is 
3 seconds for the same problera* 

Both, the saving in computation effort# and storage 
make this formulation and algorithm more attractive espeoiallLjr 
for a* large interconnected power system. The proposed 
method very easily takes the losses into account with the 
help of loss coefficients due to Kirchmayer [ 1 ,2 ] * 



GEHSHm 5 


CONCLUSION 

This chapter is aimed at reviewing the significant 
results obtained during the course of this work and making a 
few suggestions for the future line of work in this area* 

Before reviewing the work donei a brief account of the objec- 
tives of this investigation is summarised* 

As is already mentioned in chapter 1 , the problem of 
economic load scheduling of hydro*r.thermal systems is much more 
complex than that of purely thermal systems * Several attempts 
have been made in the past to solve the scheduling problem of 
hydro-thermal systems by using techniques like Dynamic programm-* 
ing [26] and NLP techniques [ 11 ,27^29^50] • They have resulted 
in solution procedures requiring large computational efforts 
and storage* With the larger interconnected power systems, 
the procedures become bogged down* 

Hence the main objective of the present work has been 
to find a suitable solution technique which can be applied to 
the economic load scheduling problem in larger interconnected 
hydro— thermal systems and which requires lesser computational 
efforts and storage than the existing methods* 

Th e proposed solution method makes use of the linearised 
objective function and constraints around a nominal point. 

In each subinterval, the successive dual LP technique is used 
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to solve the optinisation problem of hydro— thermal system* 
load domarids are taken to be deterministic* However the 
choice of water value is crucial and the efficiency of the 
method greatly depends upon it * 

The losses have been taken into account by the loss 
coefficients due to Kirchmayer[ 1 ,2 ] * The load demend 
equation (4*2) forms the equality constraint and requires to 
be strictly satisfied • Mi numerical example of short range 
scheduling problem has been taken* Based on the investigations 
made and the results obtained, the following major conclusions 
can be drawn* 

(i) The transmission losses can be taken into account 
for the load scheduling of hydro-thermal systems* 

(ii) Successive dual IP technique can be used to the 
problem of economic load scheduling* 

(iii) Sufficient reduction in the number of constraints, 
storage and computational efforts can be obtained 
by the proposed algorithm* 

BDTI3RE SOOPE OB WORK 

As a sequel to tho results obtained during the present 
investigation,; the following line of future work seems to be 
worth-pijr suing * 
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1 , For a short range problem, the Ipad demands and the dis- 
charge of water are assumed to be deterministic as the 
uncertainties in them have been ignored* But in actual 
practice, these can only be characterised probalistically 
and the scheduling problem becomes a stochastic programming 
problem* Suitable modifications of the proposed solution 
procedure may be possible to consider the stochastic nature 
of the load demands and the discharge of available water • 

2, A greater saving in the cost can be further achieved by 

optimising the reactive power of the system* Both, real end. 
reactive generations contribute to the transmission losses*. 
Hence problem should be solved by using A.0* power flow 
model for the electrical power network* A.C. power flow 
equations would represent the equality constraints of the 

electrical network and the operating limits on the voltage 
magnitudes, reactive powers and line flows etc* may be 
included in the formulation, in addition to the constraints 
of the hydro subsystems ♦ Transformer tap-settings should 
also be included in the model* 

3* The head variation of the reservoirs for a short range 

problem has been ignored in the present work. The dynamic.' 
of the reservoirs should be taken into account for the head 
variations in the optimal scheduling of hydro-thermal systomr' - 
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